Atom cooling using the dipole force of a single retroflected laser beam 
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We present a mechanism for cooling atoms by a laser beam reflected from a single mirror. The 
cooling relies on the dipole force and thus in principle applies to arbitrary refractive particles in- 
cluding atoms, molecules, or dielectric spheres. Friction and equilibrium temperatures are derived 
by an analytic perturbative approach. Finally, semiclassical Monte-Carlo simulations are performed 
to validate the analytic results. 
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I. INTRODUCTION 

Optical cooling of atoms has come a long way since it 
was first proposed; magneto-optical traps are even found 
in undergraduate laboratories [U, 0, Q . The field of ul- 
tracold molecules, in contrast, is still in its infancy. Ul- 
tracold diatomic alkali molecules (< 100 pK for Rb2 
are routinely produced from Bose-Einstein condensates 
through Feshbach resonances. Some groups (see, for ex- 
ample, and 3) have demonstrated the possibility of 
cooling the internal degrees of such molecules, cooling ul- 
tracold di-alkali molecules to their lowest rovibronic lev- 
els by means of laser-stimulated state transfer processes. 

Present methods of producing ultracold samples suffer 
from one of two major drawbacks: either they are specific 
to particular species, or they produce very low densities. 
The bulk of optical cooling methods are applicable only 
to a handful of species because they rely on a closed opti- 
cal transition within which the population can cycle 0. 
Most atoms and molecules do not have such a transi- 
tion available, but instead exhibit a large number of loss 
channels through which the population is gradually lost, 
halting the cooling. Samples of cold molecules are there- 
fore generally produced by capturing the low- velocity tail 
of the Maxwell-Boltzmann distribution of a hotter initial 
sample However, such filtering methods do not lead 
to an increase of the phase-space density and thus only 
capture a small fraction of the initial population, leading 
to very dilute samples. 

One possibility to solve this problem is through the use 
of non-resonant processes [1] or cavities [l3. El, HI, HI, 
[T^ . The latter require extremely high precision align- 
ment of the cavity mirrors as well as complicated loading 
of the molecules into the optical cavity mode. The re- 
quirements for integrated systems near the surface of a 
substrate in the form of atom chips [Tsj are even more 
stringent. 

Here, we investigate a mechanism for the cooling of 



a particle using only a single plane mirror in place of a 
cavity. In principle, this scheme only relies on the dipole 
force of a refractive particle in a laser beam and thus 
applies to a wide range of atomic and molecular species 
as well as, for example, dielectric micro- or nanospheres. 
However, in the present work we focus on the basic prin- 
ciples of the cooling scheme and thus restrict the analysis 
to the simplest case of a two-level atom. 

Conceptually, one can view the interaction between the 
atom and the mirror as being closely related to that be- 
tween a micromechanical oscillator, acting as a mobile 
mirror, and a second mirror. Such schemes have been in- 
vestigated both theoretically 0, [l3| and experimentally 
[ll,[l^ in various configurations. 

Although we have recently shown that one can treat 
these two situations as two opposite limits of the same 
model [13], the situation we explore here behaves differ- 
ently, and this can be attributed to two facts. Firstly, 
the coupling strength between the static and moving 
scatterer (atom or mirror) is very different in the two 
cases: an atom merely perturbs the field it interacts with, 
whereas a mirror acts as a moving boundary condition 
and changes the field significantly. Secondly, the effect 
we investigate here is only dominant at large atom-mirror 
separations. Thus, our proposed cooling scheme operates 
in a parameter regime that is as yet mostly unexplored. 

This paper is structured as follows. In the next section 
we introduce the key features of our system and propose 
a simple classical explanation of the cooling scheme. In 
section Hill the relevant quantum equations of motion are 
introduced. Section [IVI solves these equations of motion 
analytically through the use of perturbative methods, 
whereas section |V] gives the results of numerical simu- 
lations used to explore the implications of the theory in 
further detail. SectionlVllcompares our results with those 
of traditional Doppler cooling, and finally section I VIII 
summarizes and concludes our discussion. 



II. MOTIVATION 
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We start with a classical explanation of the situa- 
tion, which provides the motivation for the mathematical 
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FIG. 1: Schematic of the coohng scheme, (a) In the lab- 
oratory frame, an atom moving with velocity v parallel to 
a mirror, a distance rc away, interacts with a pump beam 
and its time-delayed reflection, (b) In the frame of the atom, 
the (relativistically transformed) pump beam and its reflec- 
tion are tilted by an angle 6 — arcsin(t;/c) and produce a net 
retarding force, F. 




FIG. 2: Schematic representation of the key components of 
the system under consideration. The atom is separated from 
the mirror by a distance rc and Ues in a standing wave main- 
tained by the pump beam. 

"position dependent phase locking" described in (26| . 



model presented in the next section. The phenomenon of 
optical binding has been known for some time (see, for 
example, [2ll. [2^. ) and is now a common occurrence. 
At a basic level, optical binding takes place between two 
dielectric spheres when one sphere focuses the light onto 
a second sphere, which is subsequently trapped. If we 
now consider just one such sphere in front of a mirror, 
the modified electric field will be reflected back towards 
the sphere itself. In essence, then, the sphere will be at- 
tracted to its own image. However, this interaction is 
delayed by the time 2t it takes the disturbance in the 
electric field to travel from the sphere to the mirror and 
back, where r is the time that light from the atom takes 
to reach the mirror. Suppose, now, that the sphere is 
moving parallel to the plane of the mirror with velocity 
V, as shown in Fig.[TI^a). In this case, the sphere moves a 
distance 2vt during the light roundtrip time. Thus, the 
disturbed light field lags behind the particle and creates 
an attractive force in the direction opposite to the mo- 
tion of the particle. This attractive force can be shown 
to be a viscous force, i.e., it is proportional to v. This 
scenario also applies to the case of a single atom inter- 
acting with an off-resonant beam, where the atom can 
effectively be modeled as a refractive particle. Note that 
the interaction between a single atom and its image in a 
distant mirror has already been demonstrated [2j, 25] . 

Alternatively, we may consider the same situation in 
the reference frame of the moving particle, as shown in 
Fig. [Jb) . In this case both the incident laser beam and 
its reflection are tilted by an angle 9 = arcsin(w/c) with 
respect to normal incidence on the mirror. The sum po- 
tential is therefore offset from the particle position by an 
amount proportional to its velocity, leading again to a 
velocity-dependent force opposing the motion. 

Similar arguments apply for a particle moving along 
the direction of the pump beam, i.e., orthogonal to the 
mirror. In this case, the phase of the reflected beam is 
determined by the interaction of the particle with the 
pump at the earlier time 2t. In effect, the atom exerts 
a delayed phase change on the light field, dragging the 
potential along with itself while moving and thus creating 
a non-conservative force. This effect is similar to the 



III. MATHEMATICAL MODEL 

In order to analyze the principles of the proposed cool- 
ing scheme most clearly, we simplify the situation de- 
scribed above to a one-dimensional (ID) scheme and as- 
sume a single two-level atom as the particle to be cooled. 
A schematic of the system is shown in Fig. [2] 

The atom has a transition frequency Wa and a decay 
rate 2r and is described by the operators p and x asso- 
ciated with the atomic momentum and position, respec- 
tively, and by the atomic dipole raising (ct"*") and lower- 
ing (iT~) operators. The atom is coupled to a continuum 
of quantized electro-magnetic modes with frequencies uj 
and standing- wave mode functions f{u>,x) = sin(cLia;/c), 
described by the field annihilation aiuj) and creation op- 
erators a){uj). The mirror is at position a; = 0. For 
simplicity, we neglect the frequency-dependence of the 
atom-field coupling and assume a single coupling coeffi- 
cient g. Finally, mode uj^ is pumped by a laser, which 
enters our analysis as an initial condition, and far off res- 
onant pumping is assumed, | A| = — wqI ^ F, where 
the atom mainly acts as a refractive particle and sponta- 
neous scattering is reduced. For the numerical examples 
given in this paper, we consider ®^Rb atoms and a realis- 
tic pump beam that is detuned from the 5Si/2 ~^ 5P3/2 
transition of *^Rb by several linewidths. 

The starting point for describing the coupling between 
the atom and the field modes is the quantum master 
equation: 



H,p 



Cp, 



(1) 



where p is the density operator of the full system com- 
prising all modes and the atom. Applying the dipole and 
rotating wave approximations and in a frame rotating 
with the driving frequency uj^, the Hamiltonian H reads 



2m 



h(uj — Ci;o)a^(t^)a('^) dw 



ihg / [a^ a{uj) f {uj , x) — h.c] du), 



(2) 
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and the damping term associated with atomic decay into 
modes other than the ID system modes reads 



Cp = -r 



pa a 



N{u)a 



-e-™*/5e™*a+ du 



(3) 



Here N{u) describes the ID projection of the sponta- 
neous emission pattern of the atomic dipole. In the low 
saturation regime we can adiabatically ehminate the in- 
ternal atomic dynamics and formally express the dipole 
operator as 



Y2 



9 / fi^^,x)a{uj)du; + C 



(4) 



where ^~ is a noise term (27[. 

In the following, we present two different ways of pro- 
ceeding from this point. In the first instance we approx- 
imate further and use perturbation theory to derive the 
force experienced by the atom analytically, in section ITVl 
and in the second instance we derive semiclassical equa- 
tions of motion. The latter approach is then applied to 
numerical simulations of the situation in section fVl 



We now assume coherent states at all times for the fields 
and replace the operators with their respective expecta- 
tion values. Since we are pumping the atom at a sin- 
gle frequency, we take the initial condition a{Lu,0) — 
AS{uj — Wo), where A is the amplitude of the pump field, 
such that \A\'^ is the pump power in units of photons per 
second, and S is the Dirac (5-function. We now expand 
the fields a{u!,t) in the weak-coupling limit in powers of 
the coupling constant. 



(6) 



with a„(a;, t) being the nth coefficient of the series expan- 
sion. Solving Eq. ([S]) by perturbation theory then yields 
the zeroth order term in g^D{A) 



ao{cj,t) = AS{uj - Wo), 
and the first order term 



(7) 



. ^ exp [-i{u; - ujojt] - 1 / N • / / N 

ai(uj,t) — A s\v\lox I c) swyLo^x I c) . 



We now proceed to find, to second order in g^D{A), the 
static force, F{xQ,t) = —dH/dx, acting on the atom: 



IV. ANALYZING THE MODEL: A 
PERTURBATIVE APPROACH 

A. Friction force 

We first derive an analytical approximation for the fric- 
tion on the atom. To this end, we treat atomic motion 
semiclassically and thus replace the operator x by an 
atomic position x. After inserting Eq. ^ into Eq. ([2|) 
and Eq. we obtain the Hamiltonian 

H = J — ti;o)'^^('^)£*(w ) dw -I- ng^D{A) 

X J J sin(cjia;/c) sin(ijj2x/c)a^(wi)a(ijJ2) dwi dijJ2, 

where we have defined D{A) ~ A/(A^ +r^). As a conse- 
quence of the assumed large pump detuning A, we in the 
following neglect that part of the decay term Cp which 
leads to spontaneous scattering of photons between the 
quantized modes by the atom. 

Let us first consider a stationary atom at a fixed po- 
sition X — xq. Starting with the Heisenberg equation of 
motion for the annihilation operators 

—a[uj,t) = - ^H,a{uj,t}^ , 

we arrive at the integro-differential equation 

-a(w, t) — — i{u) — Wo)a(w, t) 



dt 



ig'^D{A)sin{LJx/c) / sin(ti;ia;/c)a(ci;i, t) dwi. (5) 



F{xo,t) = - \A\^ g^D{A)ujQ{ sm{2ujoxo / c) 



-g^D{A)sm'^iujQXo/c) U cos^ {ujqXq / c) - l] 



The first term in the above equation describes the in- 
teraction of the atom with the unperturbed pump field, 
whereas the second term is the lowest order correction of 
the force due to the back-action of the atom on the light 
fields. Note that this force is independent of time. 

Similarly, we can now calculate the force on an atom 
moving at a constant velocity v. For this we assume that 
the atom follows a trajectory given by x{t) = xq + v{t — 
to), where to is a long enough time for the system to 
reach a stationary state, i.e., to is larger than twice the 
propagation time t = xq/c of the light from the atom to 
the mirror. We can then solve Eq. ([5]) up to first order in 
both ti and g^D{A). The friction force in the longitudinal 
direction is finally obtained as 

F||(xo,i) = ^^v\A\'[g^DiA)]\in'{2koXo) 



c 

2TThuJn 



^vT\Af[g^D{A)Y sm{4koxo).{8) 



The second term in Eq. ^ is larger than the first by a 
factor of the order fco^^o = ^o^o/c and is therefore domi- 
nant if the distance of the atom from the mirror is much 
larger than an optical wavelength. We may then approx- 
imate the longitudinal friction force by 



g^A 



r2 



sin(4fcoXo). (9) 
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suming a Gaussian mode function of waist w = 0.7|i.m. 
Note that F\\ and F±_ are comparable in magnitude for 
the parameters chosen here where the mode waist is com- 
parable to the optical wavelength, and that there exist 
regions where both these forces promote cooling. 

In the remainder of this paper, however, we will con- 
centrate on a one-dimensional treatment of the problem 
and therefore only consider the longitudinal friction force. 
This could correspond, for example, to the imaging ar- 
rangement of Eschner at al. p3 |. 

In terms of more familiar parameters, we can 
rewrite Eq. © in the limit |A| > T as 



FIG. 3: (Color online) Spatial dependence of the longitudinal 
heating coefficient T (thick solid line). The shaded areas pro- 
mote cooling in the longitudinal direction. Also drawn is the 
transverse heating coefficient (thin solid line) and the field 
amplitude (dotted line). Parameters are for Rb atoms and 
|Ap = 62.5r/(27r), A = -lOr, t = 0.25/r, w = 0.7 ^m. 

Supposing that the species we are cooling is ®^Rb, and 
setting |Ap = 62.5r/(27r), A = -lOF, t = 0.25/r, and 
Gaussian beam waist w = 0.7 |J^m, Eq. ^ predicts 1/e 
cooling times of the order of 2 ms. The value for r that 
we use implies a separation between the atom and the 
mirror of the order of several metres. We suggest that 
this problem can be overcome through the coupling of the 
light into an optical fibre, thereby avoiding the effects of 
diffraction. A recent experiment making use of a similar 
technique is described in [28| . 

Equation ([9]) indicates an exponential decay or increase 
in velocity. We define the heating coefficient T of an en- 
semble of atoms as the proportionality constant in the 
relation dp^ /At — Tp^, which thus depends on position 
as sin(4fcoXo). Moreover, since (x T for a thermal en- 
semble, we also have dT/dt = TT. Figure [3] shows a plot 
of T against atomic position, where we introduced the 
coordinate Xq relative to the nearest node of the stand- 
ing wave pump. It is only in certain intervals that we 
expect the longitudinal force to be a damping force, as 
indicated in this figure by the shaded regions. 

We now derive the friction force in the transverse di- 
rection, i.e., orthogonal to the pump beam. In this 
case the coupling constant g becomes a function g{r), 
where r is the coordinate in the transverse direction. 
For an atom moving at small constant velocity, we may 
write g[ro + v{t - io)] ~ giro) + v{t - tQ)g'{rQ) where 
g'{r) = dg/dr. Substituting this in Eq. ([5]) we can derive 
an expression for the friction force, F± {xq, t) ^ —dH/dr, 
in the direction of r: 

F^xoM - -47rhvT\A\' (^TT^J 

X sin^(fco2;o) cos(fcoa;o)• 
This transverse friction force is also shown in Fig. [3] as- 



F»(xo,t) = -4usr— ^fi./c^Tsin(4fco2;o), (10) 

where s = g'^\A\'^ /{A'^ + F^) is the maximum satura- 
tion parameter of the atom in the standing wave, da = 
3A^/(27r) is the atomic radiative cross-section, and where 
we used the relation 2Tr g^/T = Aaa/inw"^). 

Aside from allowing us to make predictions of cooling 
times, Eq. ^ and Eq. pH)) also highlight the dependence 
of this cooling effect on the variation of certain physical 
parameters. In particular, i^y depends on the square of 
the detuning, which means that it is possible to obtain 
cooling with both positive and negative detuning. The 
friction force also scales with w~'^ and \A\'^. Hence, for a 
fixed laser intensity, proportional to jAp/w^, i.e., fixed 
atomic saturation, friction still scales with w~'^ and thus 
a tight focus is needed in order to have a sizeable effect. 
A very promising feature of these two equations is the 
linear dependence of the cooling rate on r. In section fVl 
we further analyze the dependence of the cooling rate 
on the various parameters and support the validity of 
the analytic solution by comparing it with the results of 
simulations. 



B. Localizing the particle 

Equation ^ shows that, in order to observe any cool- 
ing effects, we need to localize the particle within around 
A/8. This can be achieved, for example, by an additional 
far-off resonant and tightly focused laser beam propagat- 
ing parallel to the mirror forming a dipole trap centered 
at a point xq. We characterize this trap by means of its 
spring constant fcj, such that the trapping force is given 
by i^t = —kt{x — xo), or equivalently by the harmonic 
oscillator frequency Ut = \/ {kt/m), where m is the mass 
of the atom. 

If we now assume that the atom oscillates as x{t) — 
Xq + sm[ujt{t — to)] in the trap with a maximum dis- 
tance Xin of the particle from the trap center and a cor- 
responding maximum velocity v^, it is possible to derive 
a new friction coefficient by perturbation theory in Xm- 
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Proceeding along the lines of section IIV Al we arrive at 
i^||(a;o,t) = — 2TrhkQVmT smc{2uJtT) 

X \A\' (^1^) sin(4fcoXo) . (11) 

Note that this formula reduces to Eq. ^ in the limit of 
small Wt- The sinusoidal dependence on uJtT can be ex- 
plained in an intuitive manner: the effect on the particle 
is unchanged if the particle undergoes an integer number 
of oscillations in the round-trip time 2t. 

While Eq. (fTTj) was derived for an oscillating particle, 
it is still only correct to lowest order in Xm and therefore 
does not include the effect of a finite spatial distribution. 
In order to obtain an estimate for the friction force in the 
presence of spatial broadening, we calculate the overall 
energy loss rate experienced by the particle in terms of 
the time average of Eq. 

/dp'\_ 2nklpl ,( ff^A Y 

X / sin[4fcoa:o + 4fcoa;niSin(r)]cos^(r)dr, (12) 
Jo 

where po — mvm is the maximum momentum of the par- 
ticle in the trap given by po = XraV^nkt- The value of 
the integral in can be expressed as 

— [sin(4A:oa;o) Ji(4/coa;m) + cos(4fcoa:o)i/i(4fcoa;m)], 

4A;oa;„i 

where Ji is the order- 1 Bessel function of the first kind 
and Hi is the order- 1 Struve function [2§|. At the point 
of maximum friction, x'q — — 3A/16, the integral in the 
above equation reduces to 27rJi(4A:oa;m)/(4fcoXm) which 
can be readily evaluated. 

For small values of Sm, the effect of this averaging pro- 
cess is to introduce a factor of one half into Eq. pT|) . 
which can be seen as being equivalent to the effect of 
cooling merely one degree of freedom when the atom is 
in a harmonic trap. Finally, Eq. (|12p is modified similarly 
to Eq. (fTTj) to include the effects of the harmonic oscilla- 
tion by replacing r sin(2ajtT)/(2wt). This results in an 
approximate expression for the friction, taking into ac- 
count the periodicity in the time delay as well as spatial 
averaging effects, 

(F||(xo,i))= - hklv^T\A\^ (a^^) ^"^^(2a.tr) 

X / sin[4fcoa;o + 4fcoa;,„sin(r)]cos2(T)dT. 

(13) 

C. Capture range 

As discussed above, the addition of the dipole trap 
introduces several features into the friction force. Plot- 
ting the variation of the friction force in Eq. (fT^ with 
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FIG. 4: Dependence of the heating rate (dp^ /dt) on the square 
of the initial momentum, pg, for ujt = 0.45 x and x'q = 
3A/16. Other parameters are as in Fig. El Cooling is achieved 
only for a finite range of initial momenta. 

the particle's initial momentum, as in Fig. [3J shows that 
the force changes sign for high enough initial momentum. 
This is due to the broader spatial distribution for faster 
particles in the harmonic trap. For fast enough veloc- 
ities, the particle oscillates into the heating regions, as 
shown in Fig. [3l even if the trap is centered at the posi- 
tion of maximum cooling. This defines a range of initial 
momenta, starting from zero, within which a particle is 
cooled by this mechanism; faster particles are heated and 
ejected from the trap. Note that this result was derived 
from the friction to lowest order in velocity and higher 
order terms are expected to affect the capture range fur- 
ther. 

At particular values of Xp, e.g. at — 3A/16, this capture 
range can be conveniently estimated by using the location 
of the first zero of the Bessel function, 

Po ~ Q.9b^^/{mkt)lkQ = 0.958muJt/kn. (14) 

Thus, Pq (X uj'^, and the capture range as defined in Fig. [4] 
is expected to scale with the square of the trap frequency. 
We compare this later in section |V] with the results of 
numerical simulations. 



D. Diffusion and steady-state temperature 

In the preceding discussion we found a friction force 
which cools an atom towards zero momentum. In prac- 
tice, the cooling process is counteracted by momentum 
diffusion due to spontaneous scattering by the atom of 
photons from the pump beam into other electromagnetic 
modes and between the two counterpropagating compo- 
nents of the standing wave pump itself. In a simplified 
Brownian motion model, this diffusion introduces a con- 
stant in the equation for dp^/dt, resulting in a constant 
upward shift of the curve in Fig. 21 This slightly re- 
duces the capture range for fast particles, but its main 
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FIG. 5; Calculated steady-state temperature Tm for an atom 
confined in a harmonic trap as a function of position whilst 
keeping the detuning and pump field constant. tJt = 0.1 x 2nr; 
other parameters are as in Fig. [S] 



effect is to introduce a specific value of the momen- 
tum where friction and diffusion exactly compensate each 
other. This point corresponds to the steady-state tem- 
perature achievable through the coohng mechanism dis- 
cussed here. 

To lowest order in the coupling coefficient g^, the dif- 
fusion is given by the interaction of the atom with the 
unperturbed, standing-wave pump field. In this limit, 
diffusion in our system is therefore identical to that of 
Doppler cooling [s^, [Ml, [H, [s^l , where the diffusion co- 
efficient D is given to lowest order in s by 



D = h^klVs 



cos^(fcoa^o) + I sin^(fcoa;o) 



(15) 



The steady-state temperature Tm of mirror-mediated 
cooling is then obtained from UbTm = Dvm/ F\\{xo,t) 
where F\^{xo,t) is the friction force given by Eq. (fTTj) . 
For lAI >> r we find 



1 h wtF 2 3cos2(fcoa;o) 
bn ks sin(2wtT) sin(4A:oa;o) 



(16) 



An example of the dependence of Tm on the trap position 
is shown in Fig. [5l predicting a minimum temperature of 
the order of 400 \xK. Whilst this may seem large in com- 
parison to the Doppler temperature of 141 (a.K, one has 
to keep in mind that Tm, given by Eq. (jl6[) . is insensi- 
tive to detuning and, for far off-resonant operation of the 
order of tens of linewidths, it will be the dominant mech- 
anism. This is further discussed in section IVTl We also 
note that Fig. [5] further highlights the importance of the 
requirement for localizing the particle. Using Eq. pU]) 
we can approximate the steady-state temperature at the 
point of maximum friction by 



T 8aa 



It is interesting to note that this expression is closely re- 
lated to the expression for the limiting temperature in 



Doppler cooling, ksT = HT, but where F is replaced 
by the inverse of the atom-mirror delay time, 1/r, and 
where a geometrical factor related to the mode area di- 
vided by the atomic cross section is included. 



V. NUMERICAL SIMULATIONS 

In this section we now investigate a more accurate 
numerical model to corroborate the simplified analyti- 
cal results obtained above. In order to render the prob- 
lem numerically tractable, the continuum of modes is re- 
placed by a discrete set of modes with frequencies ujk, 
k = 1, . . . , N . The master equation ^ is then converted 
by use of the Wigner transform into a Fokker-Planck 
equation for the atomic and field variables. Applying 
a semiclassical approximation and restricting the equa- 
tion of motion to second-order derivatives, one arrives 
at an equivalent set of stochastic differential equations 
for a single atom with momentum p and position a; in a 
discrete multimode field with mode amplitudes ak [3^ . 



dx — —dt, 



(17a) 



34| 



dp = t^[£{x)££*{x) - £*{x)££{x)]dt 

-Uo[£{x)££*{x)+£*ix)^£ix)]dt 
-kt{x-xt)dt + dP, (17b) 
dafe = iAkakdt - {iUo + j)£{x)f^{x)dt + dAk, (17c) 

where fk{x) = sin(cjfcx/c) are the mode functions, £{x) = 
X]fc o.k.fk{x) is the total electric field, = ujQ — uJk is the 
detuning of each mode from the pump, Uq is the light 
shift per photon, and 7 is the photon scattering rate 
The terms dP and dAk are correlated noise terms 
responsible for momentum and field diffusion. 

In the following, we set the trap center to Xg = 
— 3A/16, which is the point where the analytic solution 
predicts the maximum of the damping force. We use 
N = 256 field modes with a mode spacing of F/10. At 
the start of every simulation, all field modes are empty 
with the exception of the pump mode which is initialized 
at 625 photons, corresponding to a laser power of around 
50 pW for our chosen parameters. 

The simulations were performed in runs of several 
thousand trajectories. Each such run was performed at 
a well-defined initial temperature, with the starting mo- 
menta of the particles chosen from a Gaussian distribu- 
tion, and the starting position being the center of the 
trap. 



A. Friction force and capture range 

Fig. [6] presents the results of a set of simulations per- 
formed when setting the noise terms dP and dAk in equa- 
tions p7ap to zero, i.e., neglecting momentum and pho- 
ton number diffusion. The simulation data are compared 
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(a) 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.00x 10'' 
(Initial momentmn, po}^ (h'^k^) 

(b) 

FIG. 6: Comparison of heating rate (dp^/dt) for the simula- 
tions without noise (data points) with the analytic approx- 
imation, Eq. (I13|l . including the harmonic trap (solid line), 
(a) Weak harmonic trap, ut — 0.3 x 27rr, showing also the 
linear dependence in the limit of small momenta, Eq. 
(dotted line), (b) Stiff trap, ujt = 0.5 x 2nr. The trap posi- 
tion x'o = — 3A/16 and other parameters are as in Fig. O 

with the result of the perturbative calculations Eq. 
For modest values of Wt, Fig. [S^a) justifies the averag- 
ing process used to derive p2)) which was based on spa- 
tial averaging but neglecting higher order terms in v. In 
contrast, for larger trap frequencies, the numerical sim- 
ulations diverge significantly from the analytic result, as 
can be seen in Fig. ^h). We expect that the terms in 
higher powers of the initial speed, which were dropped 
in the perturbative solution, are responsible for this dis- 
crepancy. 

We have already seen, in Eq. ([14]), that the capture 
range is expected to scale as w^. For weak traps, as 
shown in Fig. [71 the numerical simulations agree well 
with these expectations. For stiffer traps, however, the 
capture range is consistently larger than that predicted; 
in fact, the simulations predict a capture range of around 
450 mK for a trap frequency of 0.5 x 27rr 



B. Steady-state temperature 

The next step in our investigation was to run simu- 
lations involving the full dynamics given by Eqs. (|17ap 
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FIG. 7: Capture range extracted from the simulations (data 
points) as compared to the analytic solution (dotted line) for 
various values of cot. The solid line is a quadratic fit to the 
data for cot > 0.3 x 2nr and is only intended as a guide to the 
eye. Other parameters are as in Fig. |6] 




-4 I ' ' ' ' ■ ' ' ' ' ' ' ' ■ ' ' ' ' ' ■ J 

0123456789 10 

Initial temperature, Tq (niK) 

FIG. 8: Heating rate (dT/dt) extracted from the simulations 
starting at a number of initial temperatures. The solid line 
represents a quadratic fit to the data, uit ~ 0.5 x 2nT; other 
parameters are as in Fig. |6] 

including the diffusion terms. Because of the discrete 
nature of the field modes with uniform frequency spac- 
ing used in the simulations, the numerically modeled 
behavior is always periodic in time with a periodicity 
given by the inverse of the frequency spacing. The sim- 
ulations therefore cannot follow each trajectory to its 
steady-state. Instead, simulations were performed in sev- 
eral groups of trajectories, each group forming a thermal 
ensemble at a well-defined initial temperature. For each 
such group of trajectories the initial value of dT/dt was 
calculated. The results for wt = 0.5 x 27rr are shown in 
Fig. [HI where the error bars are due to statistical fluctu- 
ations for a finite number of stochastic integrations. The 
steady-state temperature is that temperature at which 
dT/dt = as clearly illustrated in this figure. For the 
chosen parameters, our data suggest a steady-state tem- 
perature of 722±54 (j-K with a 1/e cooling time of around 
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Trap frequency, lj^ (27rr) 

FIG. 9; Steady-state temperature for a number of simulations 
(circles) compared to the analytic formula (|16p (solid line). 
The solid square represents the equivalent data from Fig. [8l 
resulting from a much larger number of simulations. Param- 
eters are as in Fig. [S] 
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FIG. 10: Comparison between the calculated steady-state 
temperatures for mirror-mediated cooling Tm (dash-dotted 
line), Doppler cooling Td (dashed), and in the presence of 
both effects T (solid), drawn as a function of detuning whilst 
keeping the saturation parameter constant, cjt = 0.1 x 27rr; 
other parameters are as in Fig. [G] 



3.0 ms. This compares reasonably well with the steady- 
state temperature of 597 ^iK predicted by Eq. p^ . 

We finally performed a large number of simulations 
to investigate the dependence of the steady-state tem- 
perature on the trap frequency. Equation (jl6p indicates 
that as one decreases ujt the steady state temperature de- 
creases. This is clearly seen in Fig. [HI which compares the 
prediction of Eq. (|16p with a set of numerical simulations. 
The trend in the data is reproduced well by the analytic 
expression. However, the simulated steady-state temper- 
ature is consistently a little higher than predicted. We 
expect that this discrepancy is due to one of two reasons, 
(i) Equation was derived from the friction Eq. pT|) . 
i.e., without the spatial averaging of (jl3p which would 
reduce friction, (ii) Higher order terms in the velocity 
V are also expected to reduce friction compared to the 
lowest order analytical result. In both cases, therefore, 
the analytic expression is expected to overestimate the 
friction force and thus to predict too low temperatures. 



VI. BEYOND ADIABATIC THEORY 

All the theoretical analysis and simulations discussed 
so far have been based on adiabatic elimination of the in- 
ternal atomic degrees of freedom, and therefore neglected 
Doppler cooling. In Fig. [TOl we explore the variation of 
Tm and the Doppler temperature, Td, as a function of de- 
tuning from resonance when the particle is at the point 
of greatest friction {xq = 3A/16), where To is given by 
Td = -ftr(A2 -I- r2)/(2A) for negative values of A. In 
the presence of both cooling effects, the stationary tem- 
perature achieved by the system is given by 



Thus, for the parameters of Fig. [21 the calculated steady- 
state temperature T reduces to 250 (a.K in the limit of 
vanishing cjt- 

From Fig. [10] one can see that the mirror-mediated 
force, for our tightly focused pump, is stronger than the 
Doppler force for detunings larger than around lOF in 
magnitude. In practice this has two implications: for 
large negative detunings, we expect the steady-state tem- 
perature of the system to be significantly lower than that 
predicted by Doppler cooling; whereas for large positive 
detunings, we still predict equilibrium temperatures of 
the order of mK. 

Both our pcrturbative expressions and our simulations 
are calculated to lowest orders in the atomic saturation. 
However, it is well known that in the limit of very large 
detunings also higher order terms in the saturation pa- 
rameter s become significant. Using the full expression 
for the diffusion constant [3(| , we can estimate the detun- 
ing for which we expect minimum diffusion and temper- 
ature. For the value of the saturation parameter s < 0.1 
used throughout this paper, it can be shown that Tm at- 
tains a minimum at detunings of up to several tens of 
linewidths. Our chosen parameters are therefore within 
the range of validity of the model. 

VII. CONCLUSION 

We have presented a mechanism for cooling particles 
by optical means which is based fundamentally on the 
dipole interaction of a particle with a light beam and 
therefore does not rely on spontaneous emission. The 
particle is assumed to be trapped and is simultaneously 
driven by an off-resonant laser beam. After the inter- 
action with the particle the beam is reflected back onto 
the particle by a distant mirror. The time-delay incurred 
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during the light round-trip to the mirror and back is ex- 
ploited to create a non-conservative cooling force. 

The system was analyzed using stochastic simulations 
of the semiclassical equations of motion representing a 
single two-level atom coupled to a continuum of elec- 
tromagnetic modes. The results of these computations 
were found to agree with the expectations of a perturba- 
tive analysis. Our models predict sub-mK steady-state 
temperatures for ^^Rb atoms interacting with a tightly 
focused laser beam several meters from the mirror, in 
an arrangement similar to that of Ref. [Ill . While most 
of the theory is presented for a one-dimensional model, 
results for the friction force in the transverse direction 
suggest that three-dimensional cooling is possible with 
this scheme. 

The model presented here requires a large separation 
between the atom and the mirror, of the order of several 
meters, for an observable cooling effect. This limitation 
can be overcome in several ways. First, the light could 
propagate in an optical fiber between the atom and the 



mirror to avoid the effects of diffraction. Second, the 
required delayed reflection could be achieved through the 
use of a cavity instead of a mirror; in contrast to cavity- 
mcdiatcd cooling schemes [l^, [HI [H, [3 13 , the atom 
would remain external to the cavity. For a time delay t 
of order 1 ns one would require a cavity quality factor 
Q = LOT [3^ of the order of 10^ — 10^, which is achievable 
with present-day technology [36i |. 
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